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Abstract 

In this paper we revisit the results of Loynes (1962) on stability of queues for ergodic 
arrivals and services, and show examples when the arrivals are bounded and ergodic, 
the service rate is constant, and under stability the limit distribution has larger than 
exponential tail. 
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1 Introduction 



The analysis of a queueing model consists of two steps: stability study and the characteriza- 
tion of the limit distribution. In this paper we consider long range dependence, i.e. assume 
ergodic arrivals and services. Concerning stability revisit the result of Loynes (1962), who 
extended the Markovian approach in an elegant way. For the weak dependent situation the 
limit distribution usually has an almost exponential tail. Here we show counterexample for 
ergodic situation, i.e. when the arrivals are bounded and ergodic, the service rate is constant. 

2 Stability 

2.1 General result 

Let X be arbitrary random variable. Define 

X n +i = (X n + Z n+ i) + (1) 

for n > 0, where {Zi} is a sequence of random variables. We are ineterested in the stability 
of {Xi}, i.e. we are looking for conditions on {Zi} under which has a limit distribution. 
For the classical Markovian approach {Zi} are independent and identically distributed, when 
stability of {Xi} means that there exists a unique limit distribution of X n . Here consider 
the case when {Zi} is only stationary and ergodic. 

Following Lindvall (1992) introduce a stronger concept of stability: 

Definition 1 We say that the sequence {Xi} is coupled with the sequence {X[} if {X[} is 
stationary and ergodic and there is an almost surely finite random variable r such that 

X' n = X n 

for n > t. 
Put 

Vo = 0, 

n-l 

V n = ]TZ_ 4 , (n> 1). 

i=0 

Theorem 1 If {Zi} is stationary and ergodic, and E{Zj} < ; then {Xi} is coupled with a 
stationary and ergodic {X^} such that 

X' = sup V n . 

n>0 

Proof. 

Step 1. Let X_ N _ N = and define X_ N , n for n > —N by the following recursion, 

^-JV.n+l = {X- N , n + Z n+l ) + for U > -N. 
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We show that X_at i0 is monoton increasing in N, and almost surely, 

lim X_at = X' 

where 

X' = sup V n , 

n>0 

and X' is finite a.s. 

Notice that X- NjTl+1 = (X-N, n + for n > —N. First we proove that for n > —N, 

X-N-n = max{0, 

Zni Z n -\- Z n —\ ) . . . , Z n + . . . + Z_ N+ i}. (2) 

For n = -N + 1, 

X-N-N+l = (X-N^-N + Z^N+l) + 

= (Z- N+1 ) + 

= max{0, Z_ N+1 }. 



For n = -N + 2, 

X-N-N+2 = (X^N-N+l + Z_ N+2 ) + 

= max{0, X- N - N+ i + Z_ N+2 } 

= max{0, max(0, Z_ N+l ) + Z_ N+2 } 

= max{0, Z^ N+2 , Z- N+2 + Z_ N+1 }. 

Now we proove by induction from n to n + 1. 

-X"-JV,ra+l = (^-AT,n + ^n+l) + 

= max{0, max{0, Z n , Z n + Z n _ u ...,Z n + ... + Z_ N+1 } + Z n+1 } 
= max{0, Z n+ i, Z n+ i + Z„, . . . , Z n+1 + . . . + Z_ N+ i}. 

We have completed the proof of (2). Thus 

X„ A r j0 = max{0, Z , Z + Z_ u ...,Z + ... + Z_ N+1 }, 

which imlies that X_jv,o is monoton increasing, since the maximum is taken over larger and 
larger set. It remains to prove that X_jv )0 converges to a random variable X' which is finite 
a.s.. Now by the Birkoff strong law of large numbers for ergodic sequences, a.s. 

jfe,^ £ Z t = E Zl <0 

JV i=-N+l 

(cf. Theorem 3.5.7. in Stout [10]), hence a.s. 

o 

lim Zi — — oo. 
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We got that there is a random variable r such that for alH > 

oo > X_ Tj o = X_ r _j i0 , 

and therefore 

X' = sup V n . 

n>0 



X>=X' 



Step 2. Put 

and for n > 0, 

-^n+i = (X' n + Z n+ i) + . 

We show that {X!-} is stationary and ergodic. 
For any sequence = (..., z-i, z , z±, . . .) put 

F ( z -oo) = lim max{0, z ,z + Z- 1 ,...,z + ... + z_ N }. 

Then by Step 1 

X' Q = X' = F{Zf O0 ). 
We prove by induction that for n > 0, 

X' n = F(T n Z°°J, 

where T is the left shift. For n — \, 

F{TZ^) = lim max{0, Z u Z 1 + Z , . . . , Z ± + Z + . . . , Z_ N+1 } 
= ( lim max{0, maxiO, Z , . . . , Z + . . . , Z_ N+2 } + ZA) 
= (max{0, [ lim max{0, Z , . . . , Z Q + . . . , Z_ N+2 }\ + Z x }) 

N^oo 

= K 

Now we prove from nton + 1. 

X' n+ i = (X' n + Z n+ i) + 

= (F(T n Z^J + Z n+1 ) + 

= (lim max{0, Z n , Z n + Z n _ ± , ...,Z n + Z n _ x + Z n _ N } + Z n+1 ) 

= lim max{0, Z n+ ±, Z n+ i + Z n , . . . , Z n+ i + Z n + . . . , Z n+ i_N} 

= F(T n+1 Z™ 00 ). 

Step 3. Similarly to the proof of Step 1, 

X n = max{0, Z n , Z n + Z n _i, . . . , Z n + . . . + Zi, Z n + . . . + Z\ + X }, 



+ 
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and 

X' n = max{0, Z n , Z n + Z n _ u . . . , Z n + . . . + Z 1 , Z n + . . . + Z 1 + X' Q }. 
But for large n, both 

Z n + . . . + Z 1 + X < 

and 

Z n + . . . + Z 1 + X' Q < 0, 

and so 

X n = X' n = max{0, Z n , Z n + Z n _i, . . . , Z n + . . . + Z t }. 
The proof of Theorem 1 is complete. 

Remark 1. From the proof it is clear that the following extension is straightforward: If 
{Zi} is coupled with a stationary and ergodic {Z-}, and E{Z t '} < 0, then {Xi} is coupled 
with a stationary and ergodic {X^} such that 

X' = supK, 

n>0 

where 

V = 0, 

n-l 

V n = J]^( n >l). 

i=0 

2.2 Queue length for discrete time queueing 

As an application of Theorem 1 consider a discrete time queueing with constant service rate 
s, and denote by Y n the number of arrivals in time slot n. Let the initial length of the queue 
Qo be arbitrary non-negative integer valued random variable. Then 

Qn+l = {Qn — S + Y n+ i) + 

for n > 0. Put 

V = 0, 

n-l 
i=0 

Corollary 1 // {Fj} is stationary and ergodic, and E{Fj} < s 7 ^/ien i/ien {Qi} is coupled 
with a stationary and ergodic {Q[} such that 

Q' Q = sup V n . 

n>0 
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Proof. Apply Theorem 1 for 

Remark 2. This result together with Remark 1 has some consequences for network of 
servers (tandem of queues), when the output of a server (Q n + Y n+ i — Q n +i) is the input 
of another server. It is easy to show that the output (Q n + Y n+ i — Q n+ \) is coupled with a 
stationary and ergodic sequence, and the expectations of the input and the output are equal, 
so the stability condition holds for the next server, too. 

2.3 Wating time for generalized G/G/l 

This is another application of Theorem 1. According to Lindley (1952) consider the extension 
of the G/G/l model. Let W n be the waiting time of the n-th arrival, S n be the service time 
of the n-th arrival, and T n+ i be the inter arrival time between the (n+l)-th and n-th arrivals. 
Let Wo be an arbitrary random variable. Then 

W n+1 = (W n -T n+1 + S n ) + 

for n > 0. Put 

V = 0, 

n-1 
i=0 

Corollary 2 (Extension of Loynes (1964)) If{Si-i— T{\ is stationary and ergodic, E{Si_i} < 
E{Tj} 7 then {Wi} is coupled with a stationary and ergodic {W[\ such that 

W^ = supK- 

n>0 

Proof. Apply Theorem 1 for 

7 — Q — T 

u n—l - L n- 



3 Limit distribution 

In a queueing problem the properties of the limit distribution are of great importance. In this 
section we consider the special case of Section 2.2., when the arrivals {Y n } are ergodic and the 
service rate is constant s. If {Y n } are weakly dependent then the tail of the limit distribution 
is almost exponential, which may result in efficient algorithms for call admission control (cf. 
Duffield, Lewis, O'Connel, Russel Toomey (1995)). The exponential tail distribution can be 
derived using large deviation technique (cf. Glynn, Whitt (1994)). The basic tool in this 
respect is the cummulant moment generating function: 

A(0) = lim -logE{e e £Li y fc}, 



5 



assuming that this limit exists. If the set 



{9;\(9)-9s < 0}, 

is not empty then put 

5 = sup{#; A(0) - 9s < 0}. 

Then for large q 

P{Q > q} ~ e- Sq . 

The question is whether under the stability condition E(Yi) < s one has exponential tail 
distribution. 

Glynn, Whitt (1994) gave a positive answer under weakly dependent {Y n }. The limit in 
the definition of X(9) exists only under some conditions, for example, if {Y n } form a binary 
Markov chain then X(9) can be calculated (Dembo, Zeitouni (1992)), an explicite bound on 
P{Q > q} can be given (Duflield (1994)). 

For a possible extension the other problem is that E(Yi) < s is only a necessary condition 
that the set {9; X(9) — 9s < 0} is not empty, but not sufficient. This can be seen by Jensen's 
inequality: 

ilogE{e e £Li*H > -loge eE ^i yfc > = 9E{Y 1 }, 
n n 

therefore if the set {9; X(9) — 9s < 0} is not empty then E(Yi) < s. 

For long range dependent arrivals Duflield, O'Connel (1995) proved that the tail may not 

be exponential. They introduced the scaled cummulant moment generating function: 

\*(9) = lim ^-logE{e e ^ E ^ Yfc - ns] }, 

assuming that this limit exists, where a(t) and v(t) are monoton increasing functions. If the 

set 

{9;\*(9)<0}, 

is not empty then put 

5 = sup{9;\*(9) < 0}. 

Then for large q 

P{Q>q}~ e - Sv{a ~ 1{q)) . 

Duflield, O'Connel (1995) applied this result when {Y% =1 Y k } is a Gaussian process with 
stationary increments, or Ornstein-Uhlenbeck process, or a squared Bessel process. These 
examples can be motivated by multiplexing many sources. For all these examples Y n is 
unbounded. In this section we consider bounded Y n , especially binary valued Y n . Show 
examples such that Q has larger tail than exponential. 

Proposition 1 There is a stationary, ergodic and binary valued {Y n } such that E{Yi} < 1/2 
and with 3 = 3/4 the queue length sequence is stable, and for each 5 > there is q such that 

P{Q >q}> e' Sq . 
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Proof. From Theorem 1 

Q = sup V n , 

n>0 

therefore for any n 



n-l 

P{Q >q}> P{V n >q} = P{£ Y-j >ns + q}. 

3=0 

We show that for 6 t = 7, ft = 2i 2 , m = T' 1 (i > 16), 

n;-l o 

P(E K -i>7 n * + ft)> 2 "* W - 

3=0 4 

Step 1. We present first a dynamical system given in Gyorfi, Morvai, Yakowitz (1998). We 
will define a transformation T on the unit interval. Consider the binary expansion rf 3 of 
each real-number r e [0, 1), that is, r = T{2r l . When there are two expansions, use the 
representation which contains finitely many l's. Now let 

r(r) = min{i > : = 1}. (3) 

Notice that, aside from the exceptional set {0}, which has Lebesgue measure zero r is finite 
and well-defined on the closed unit interval. The transformation is defined by 



(Tr). = 



1 if < i < r(r) 

if i = r(r) (4) 
Ti if i > r(r). 



Step 2. We show that the transformation T is ergodic. 

Notice that in fact, Tr = r - 2~ r M + E[=i )_1 2 _i ! . All iterations T fc of T for -oo < fc < oo 
are well defined and invertible with the exeption of the set of dyadic rationals which has 
Lebesgue measure zero. In the future we will neglect this set. Transformation T could be 
defined recursively as 

f r - 0.5 if 0.5 < r < 1 
1 ±±^ if0<r<0.5. 



Tr 
Let 
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&i — {-^0? • • • ? -^2 l -ll 

be a partition of [0, 1) where for each integer j in the range < j < 2* 7j is defined 
as the set of numbers r = Y^=i r v^~ v whose binary expansion 0.ri,r2, . . . starts with the 
bit sequence ji, j 2 , ■ ■ ■ ,ji that is reversing the binary expansion jj, . . . , J2, ji of the number 
j = I]; =1 2 i ~ 1 j i . Observe that in there are T left-semiclosed intervals and each interval 
Ij has length (Lebesgue measure) 2~\ Now Jj is mapped linearly, under T onto 7j_ 1 for 
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j — 1, . . . , 2 l — 1. To confirm this, observe that for j — 1, . . . , 2 l — 1, if r G i] then 



r(r)-l oo 
J=l J=T(r)+l 

= r-±2-\j l -{j-l) l ) 
i=i 

i oo 

= EO' " I)i2-' + E ri2~ l . 
i=i i=i+i 

Now if < r G Jq then r ( r ) > * an d so Tr G i^-i- Furthermore, if r G i^-i then 
n = . . . = r< = 1, and thus conclude that (r _1 r)i = . . . = (T _1 r)i = 0, that is, T~V G /j. 
Let r G [0, 1) and n > 1 be arbitrary. Then r G J" for some < j < 2™ — 1. For all 

J - (2" - 1) < k < j, 

n oo 

T k r = E(J - k\2- 1 + E (5) 

1=1 l=n+l 

Now since = Ij + i for i > 1, j = 0, . . . , 2 l — 2, and the union over i and j of these sets 

generate the Borel a-algebra, we conclude that T is measurable. Similar reasoning shows 
that T _1 is also measurable. The dynamical system J 7 , fx, T) is identified with f2 = [0, 1) 
and T the Borel cr-algebra on [0, 1), T being the transformation developed above. Take \i 
to be Lebesgue measure on the unit interval. Since transformation T is measure-preserving 
on each set in the collection {Jj : 1 < j < 2* — 1, 1 < i < oo} and these intervals generate 
the Borel a-algebra J 7 , T is a stationary transformation. Now we prove that transformation 
T is ergodic as well. Assume TA = A. If r G A then T l r G A for — oo < I < oo. Let 
i? n : [0,1) — > {0,1} be the function i? n (r) = r n . If r is chosen uniformly on [0,1) then 
i2i,i?2, ••■ is a series if i.i.d. random variables. Let JF„ = a(Rn, R n +i, . . .). By (5) it is 
immediate that A G fl^L 1 JF n and so A is a tail event. By Kolmogorov's zero one law n(A) is 
either zero or one. Hence T is ergodic. 

Step 3. We define a partition of [0,1) in the following way. Let A — 0, B — [0,2~ 2 ), 
Co = A |J B . In general, for % > 1 let 

A = 2 (J V^O^- 2 - 1 ) = 2 (J 1 If +1 (6) 
j=0 i=o 



and 



We show that 



2 ! -l 2 ! -l 



(J T-^"[0,2- ai - 1 )= (J /f +1 . (7) 

j=2 i ~ 1 3=^ x 



Since [O^" 2 *" 1 ) = J 2i+1 it is clear for < j < 2 l the sets T~i[0, 2~ 2i ~ l ) are disjoint. Thus 

f i(A i )=»(B i ) = 2-*- 1 2 i - 1 . 
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Step 4. Put 

C i = A i ]jB i . (8) 

and 

oo 

c = IJ Q. 

3=0 

We show that 



-2i-l\ 



Since A = and for % > 1, 

3=0 

-l-l _! \ / 2 >-i_ 1 

(J T-^0,2- 24 - 1 ) |J |J TJ'fO^- 2 *- 1 ) 

i=0 / \j=2*-i-i 

V 3=0 J \j=2'-i-i 

= Cj-i 

and so C Cj_i, that is, 

oo 

c = U Si. 

i=0 

Furthermore 

oo 

/i(C) < J>(A) 

i=0 

oo 

= ^2~ 2i - 1 2 i - 1 

i=0 

oo 

= E2" 1 - 2 

i=0 
1 

2' 

Step 5. Define the binary time series {Yj} as 



1 if e C 
otherwise. 



Clearly {Yj} is stationary and ergodic since the dynamical system itself was so. EYi < 0.5 
since fi(C) < 0.5. If u e A> then by (6), (7) and (8) 



LO 
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that is 

y H = i,...,y- ( ^-i)H = i. 



Furthermore for i > 16, ^ = > ^ = and so = |rij + |rij > |nj + Thus 

P(E y -^>7^ + *)>MA). 

i=o 4 

By Step 3, for i > 16, 

//(A) = 2 _i ~ 2 > 2~ 2i = 2~ (1/i)2i2 = 2 - * i9i . 
The proof of Proposition 1 is complete. 
One could define 

A(0) := limsup-logEe 9 SLi^. 

n— >oo 72 

The next poposition shows that for the stationary and ergodic time-series just defined, 

{6;\{0) -6s < 0} 

is the emty set when s — 1. 

Proposition 2 For {Yj} defined in the proof of Proposition 1, 

1 Y^ n !~ 1 flv 

lim — logEe 2 -^ 3 = 9. 



Proof. By Step 3 and 5 of the proof of Proposition 1 

1 S^"'^ 1 9Y- 

6 > lim sup — log Ee^^ 3 

i— »oo Tli 

> l im in f _J l og2 E2^"=o 1 (log2 e)f ^ 

j->oo nj log 2 e 

> lim inf— ^ log 2 (2 (log2e)en >M i )) 

j^oo m log 2 e 

= lim inf 6 H — log 2 fxiAA 

i-^oo m log 2 e 

o-i+l 

= 9+ lim- log2-^ 2 = 0. 

i^co log 2 e 

The proof of Proposition 2 is complete. 
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